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Abstract 



Based on density-functional theory we analyze the full phase diagram, 
the occurrence of long-ranged orientational order, and the structural prop- 
erties of dipolar fluids. As a model system we consider the Stockmayer fluid 
that consists of spherical particles interacting via a Lennard-Jones potential 
plus dipolar forces. For sufficiently strong dipole moments one finds a region 
where a fiuid phase with long-ranged orientational order is stable. For all 
sample shapes with the exception of a long thin needle this phase exhibits a 
spatially inhomogeneous magnetization which depends on the actual shape. 
We determine the details of the magnetization structure in a cubic sample in 
the absence and in the presence of an external magnetic field. One obtains 
a vortexlike structure with an escape of the magnetization into the axis di- 
rection near the vortex axis and two point defects where the absolute value 
of the magnetization is strongly reduced. If the spherical cores of the par- 
ticles are replaced by elongated or oblate shapes a nematic phase without 
spontaneous magnetization is also possible due to the anisotropic steric inter- 
actions. We study the interplay of this nematic ordering with ferromagnetism 
in fluids of dipolar hard ellipsoids. Orientational order arises locally in the 
isotropic fluid phases near the liquid-gas interface of the Stockmayer fiuid. 
Density-functional theory allows us to determine density and orientational 
order profiles as well as the surface tension of this interface. 
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I. INTRODUCTION 



Spontaneous symmetry breaking induced by phase transitions is one of the central themes 
of statistical physics. From a homogeneous and isotropic liquid phase transitions can occur 
to phases with positional or orientational long-ranged order of the individual particles. Here 
we are mainly interested in the case of pure orientational order. Well known examples of 
such systems are the nematic liquid crystals. Initially, it was suspected that dipolar inter- 
actions are responsible for the formation of the nematic state Later, it became clear 
that the short-ranged steric interactions between strongly anisotropic (elongated or disk- 
like) nonpolar molecules are sufficient to generate a nematic phase. On the other hand, it is 
an interesting question whether dipolar interactions alone, i.e., in the absence of anisotropic 
steric interactions, are also capable of stabilizing an orient at ionally ordered liquid. That this 
is indeed possible has been shown for the first time some years ago in computer simulations 
of a dipolar soft sphere fiuid by Wei and Patey [0,^. Similar observations were made for 
dipolar hard spheres and the Stockmayer fiuid [0, in which the spherically symmetric 
Lennard- Jones potential is added to the dipolar interactions. Depending on whether one 
deals with electric or magnetic dipole moments this new phase, in which the dipoles are 
aligned over macroscopic distances, is called a ferroelectric or ferromagnetic liquid. In the 
following we will adopt the magnetic language keeping in mind that all results are equally 
valid in the electric case. Theoretically this phase transition was predicted within mean- 
field theory for a dipolar lattice gas [^] and by a generalized van der Waals theory 0. A 
density-functional theory using as input direct correlation functions of the isotropic liquid 
obtained from an integral equation theory qualitatively confirmed these findings |T^. How- 
ever, more recent work along the same lines predicts that the ferromagnetic liquid is 



metastable with respect to freezing. It might be possible to improve upon this approach by 
directly calculating the direct correlation function in the ordered phase using the methods 
of Ref. |jT2|. The present authors also developed a density-functional description |[19|J2T[ 



us- 



ing a relatively crude approximation for the pair correlation function, whose results will be 
reviewed in this paper. In contrast to Ref. |10] and to the simulations this approach enables 
one to map out whole phase diagrams which turn out to exhibit an interesting topology 
varying with the dipole strength. Other authors have focused on the isotropic-ferromagnetic 
spinodal which corresponds to a divergence of the dielectric constant, whose position can 
be estimated using correlation functions from different integral equation theories |]13|,|14|. Fi- 
nally, we mention perturbation theory as an alternative approach to determine successfully 
liquid- vapor coexistence of Stockmayer fluids with rigid or polarizable dipoles [p!5|-p!8|] . 

In all theoretical treatments particular care must be taken of the long range of the dipolar 
interaction |T^,|2D| which may give rise to diverging integrals. Several approaches to handle 
these problems have been suggested ||T0|Jl9|j20[| , the most consistent one being an explicit 



performance of the thermodynamic limit starting from a large but flnite sample of a given 
shape ||19|. For a homogeneously magnetized fluid one flnds an energy contribution that de- 
pends on the sample shape. However, this shape dependence disappears if inhomogeneously 
magnetized conflgurations are taken into account p2|,pT|. 



These complications do not arise in Heisenberg fluids governed by short-ranged exchange 



interactions for which similar phase transitions and phase diagrams have been found |23-25 



For this model as well as for the corresponding Ising fluid ||2^ there is an ongoing discussion 
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concerning the values of the critical exponents describing the isotropic-ferromagnetic transi- 



tion p6|-pq| . A Potts fluid with q > 3 states exhibits only first-order isotropic-ferromagnetic 
transitions For a review of the properties of dipolar and spin fluids as obtained from 
simulations see Ref. . 

Experimentally a ferromagnetic liquid phase in dipolar fluids has not yet been found. 
The dipole moments of molecular fluids seem to be too low, though only by a small factor, 
compared to those that led to ferromagnetism in the simulations. In ferrofluids [^, i.e., 
colloidal suspensions of small single-domain ferromagnetic particles, chemical tailoring fa- 
cilitates higher reduced dipole moments; the relevant dimensionless quantity is fcsTm^/o"^ 
where m is the dipole moment and a the particle diameter. However, in the simulations the 
phase transition always occurs at high packing fractions that are presently not attainable in 
ferrofluids. Due to the combination of fluidity and strong paramagnetism ferrofluids already 
have a wide range of technical applications that would certainly be enlarged if a ferro- 
magnetic phase could be established. Electro-rheological and magneto-rheological fluids ||^ 
also have high technological relevance. In these systems dipole moments are induced in the 
colloidal particles by external electric or magnetic fields and large differences between the 
permitivities of the solute and the solvent which then give rise to strong dipolar particle in- 
teractions having a pronounced effect on their rheological properties. However, spontaneous 
orientational order in the absence of external fields is not possible in these fluids. 

Ferromagnetic metals become paramagnetic already far below their melting tempera- 
ture. Therefore it came as a surprise when recently ferromagnetism was detected for the 
first time in a liquid metal. In Co8oPd2o the melting temperature is sufficiently close to the 
Curie temperature so that ferromagnetism in the liquid state can be reached by undercool- 
ing [32-34]. For the formation of the ferromagnetic order in these materials short-ranged 
exchange interactions play an important role. However, in addition dipolar interactions are 
also present and, as in solid ferromagnets, are essential in forming the domain structure. 
Therefore for these systems we expect strong similarities to the behavior of dipolar fluids as 
discussed here. 

Two other closely related peculiarities of simple dipolar fluids have attracted much at- 
tention recently. Strongly dipolar hard or soft spheres form chains of head-to-tail aligned 
particles at very low densities |]35| , |36[| . The chain formation hinders the condensation of liquid 
droplets. Thus a certain amount of spherically symmetric dispersion interaction is necessary 
in order to obtain liquid-gas coexistence in a dipolar fluid ||3^ (but see also Ref. ^^). Al- 
ternatively, liquid-gas coexistence of hard dipolar particles is promoted if they are slightly 
elongated, which renders the energies of the head-to-tail and the antiparallel side-to-side 
configuration more equal. It has been shown that for a narrow regime of aspect ratios a 
liquid phase condensates The relation between chaining and condensation has been 



studied in a number of theoretical papers pZ)|-fl^, 

The solid phases of dipolar fluids are less well studied. In simulations a body-centered 
tetragonal structure, which has also been found experimentally in electrorheological fluids 
P3|| , as well as fee structures with different types of orientational order were examined 
0, but it was not possible to sort out the thermodynamically stable phase. Two early 
density-functional approaches p^ , p5| did not flnd a ferromagnetic solid, but both suffer 
from algebraic errors E3] and an incorrect treatment of the long-ranged dipolar interaction 



. Our own density-functional results for the freezing of the Stockmayer fluid are presented 
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in Sec. 0. Ferromagnetism in amorphous dipolar solids is discussed in Refs. [^H]-^. 

The ferromagnetic liquid should be more stable in fluids of oblate instead of spherical 
particles. Indeed a ferromagnetic phase is formed in a fluid of cut spheres while at lower 
temperatures a columnar phase with ferromagnetic order only within the columns is ob- 
served . The relations between nematic and ferromagnetic order in fluids of dipolar hard 
ellipsoids are discussed in Sec. within the framework of a simple density-functional theory. 

Once the bulk phases and their coexistence are established one can study the interfacial 
properties which arise if different phases are brought into spatial contact with each other 
or with external walls. These cases can, in principle, be analyzed using the same density 
functionals as for the bulk. This is demonstrated in Sec. |V11| for the interface between the 
isotropic gas and the isotropic liquid which exhibits local orientational order. For dipolar 
hard spheres and Stockmayer fluids other theoretical approaches to this problem can be 
found in Refs. ||5^-|54|. There is also an extensive literature concerning more realistic models 

JSH], integral 



of specific substances such as water or methanol using lattice-gas theory 
equation theory [^, density-functional theory ||58|, and Monte Carlo simulation [ ^ . 

Non-planar interfaces arise in liquid drops and clusters. Metastable spherical liquid drops 
have been studied with density-functional theory by Talanquer and Oxtoby in con- 

nection with the analysis of liquid nucleation in supersaturated dipolar gases. In computer 
simulations of small dipolar clusters |62,63] interestingly at low temperatures one finds vor- 
texlike orientational order similar to the bulk structure of the ferromagnetic phase as dis- 



cussed in Sec. Ill 



II. MODEL AND THEORETICAL APPROACH 



We study the Stockmayer model of a dipolar fluid. It is defined by the interaction po- 
tential 

U7(ri2, uj, u') = WLjiru) + Wdip(ri2, uj, u') (1) 

where ri2 = r — r' is the interparticle vector and the orientation of the dipole moments 
is denoted hj u = {9,(p). The Lennard- Jones potential wijir) = 4e[(a/r)-'^^ — (cr/r)^] with 
parameters e and a describes the spherically symmetric dispersion interaction as well as the 
short-ranged repulsion. The dipolar potential due to point dipoles embedded in the center 
of the particles is given by 

2 

Tfl 

Wdip{r, u, u') = -— (3(m(u;)r)(m(u;')r) - m{uj)m{uj')) (2) 

where carets denote unit vectors and m is the absolute value of the dipole moment. A 
natural extension of this model are mixtures of two polar or polarizable species which also 
are important in a number of technicological applications p4| , |65| . 

For a start the Stockmayer fluid is treated within the framework of a relatively simple 
density-functional theory, which can be derived from an exact expression for the density 
functional using the low-density limit exp(— for the pair distribution function |^6[. (For 
more sophisticated choices see Ref. []TT|.) The resulting functional of the one-particle number 
density p(r, u) of particles at point r and with orientation u is 
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n[{piv,u)},T,fx] = d'rdujp{v,uj) [In {AnA' piv,uj)) - l] + j £r fcs{p{r)) 

-^j d'rd'r'dudu' p{r,u)p{r',u')Qiru - a)f{r,2,uj,uj') (3) 
— p d^rdiup{r, lu). 



Here A denotes the thermal de Broghe wavelength, p(r) = J dujp{r, u) is the orientationally 
averaged density profile, / = exp(— — 1 is the Mayer function, P = l/{kBT) corresponds 
to the inverse temperature, and p is the chemical potential. The first term in Eq. corre- 
sponds to the ideal gas free energy while the second term represents the excess free energy 
of the hard sphere reference system in a local density approximation with the Carnahan- 
Starling free energy density of hard spheres of packing fraction r] = ^pd^: 

^fcsip) = P^^- (4) 

The temperature dependent diameter d = dr{l — exp(— /?Wij(r))), as given by Barker 
and Henderson |^^, accounts for the actually soft repulsion. The third term in Eq. 
describes the contributions to the free energy due to the long-ranged part of the interactions 
in a mean-fieldlike manner. The equilibrium density distribution for a given thermodynamic 
state minimizes the above density functional. 

For bulk fluid phases the orientationally averaged density p is taken to be constant and 
the orientational distribution to be axisymmetric, so that with 



p{r,uj) = p{r)a{r,uj), j du a{r,uj) = 1, (5) 
in general, here one has 



p(r,a;) = pa(cos6') = ^ ^ a;P,(cos 6*). (6) 



oo 



27r 



1= 



In the isotropic phase q;(cos6') = l/(47r), i.e., ai = Sio/2, while in the ferromagnetic phase 
also the higher expansion coefficients ai>i take on non-zero values. If one inserts the expan- 
sion in Eq. (|^) into the density functional and employs a suitable expansion of the Mayer 
function into rotational invariants of the three spatial angles u), u' and uJi2, one encounters a 
conditionally convergent integral that is proportional to af [|19|. This is a consequence of the 
slow decay of the dipolar potential (and therefore of the Mayer function). The value of this 
integral for fluid volume ^ oo depends on the shape of the sample adopted for this lim- 
iting process. This integral was determined for ellipsoidal shapes of arbitrary aspect ratio k 
and shown to be related to the demagnetization factor known from classical magnetostatics. 
On the other hand it has been shown by general arguments that the free energies and the 



phase diagrams of dipolar systems in zero fleld do not depend on the sample shape |^ . This 
apparent contradiction is reconciled by the observation that the true equilibrium state of the 
ferromagnetic liquid corresponds to such a shape-dependent inhomogeneous conflguration 
that in the resulting free energy the shape dependence drops out. Using a general spatially 



5 



inhomogeneous orientational distribution, within the present theory we have confirmed that 
the equihbrium free energy indeed does not depend on the sample shape and is the same 
as that of the homogeneous state in a needle-shaped volume {k oo) [^. Therefore the 



phase diagram determined for this special case is valid for all shapes keeping in mind that 
for a general shape the ferromagnetic phase exhibits a nontrivial spatial structure. Details 
of this structure are discussed in the next section. We remark that these inhomogeneities 
are surpressed by the so-called "tin-foil" boundary conditions that are usually applied in 
computer simulations, i.e., by surrounding the sample with an infinitely permeable material 
(a perfect conductor in the electric case, respectively). The free "magnetic charges" in the 
surrounding material cancel the demagnetization field which leads to shape independence 
also for homogeneous magnetization. In the following we only consider the more realistic 
open boundary conditions. 

The properties of the Stockmayer fluid can be described in terms of the following dimen- 
sionless quantities: 

.3 



T* = ksT/e, p* = pa\ m* = m/Va^e. (7) 

Figure displays the phase diagrams in the temperature-density plane as obtained from the 
present theory for three different values of the reduced dipole moment m*. At low dipole 
moments there is a large coexistence region of the isotropic gas and liquid phases and a 
second-order transition to the ferromagnetic liquid along a line of critical points. This line 
ends at the coexistence curve in a critical end point. Upon increasing the dipole moment 
the isotropic-ferromagnetic transition shifts to lower densities and higher temperatures and 
a tricritical point emerges below which the transition becomes first order. Moreover all 
three fluid phases coexist at a triple point. At even higher dipole moments the gas-liquid 
coexistence becomes metastable and is preempted by the isotropic gas - ferromagnetic liquid 
transition, so that one is left with only two fluid phases. A similar series of phase diagrams 
has been obtained from a phenomenological theory 0, while in simulations to date neither 
the order of the transition nor the topology of the phase diagram have been mapped out. 



III. INHOMOGENEOUS MAGNETIZATION 

The arguments leading to the proof of the shape independence mentioned above provide 
also a characterization of the actual inhomogeneous magnetization configuration for a given 
sample shape. Locally it should exhibit the same orientational distribution that has been 
determined for the needle shaped sample, but with a spatially varying preferential direction. 
This variation should correspond to a magnetization field M(r) which on the scale of the 



system size has zero divergence and a vanishing normal component at the surface pi 
However, for a given shape these conditions do not determine M(r) uniquely and we are 
not aware of a systematic method to construct fields consistent with them. Therefore we 
pursued an alternative approach by computing this field via a quasi-free minimization of the 
density functional. A cubic discretization lattice with lattice constant a and lattice vectors 
R is superimposed on a cubic sample which we have chosen as an example. From Eq. one 
derives the following approximation for the free energy difference AF between the isotropic 
and ferromagnetic phases ||68| | : 
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3AF = pa' 5: E (i'f^Znl'+l) + l"'"' ^ E M:(R)u.«(R - R')A/,(R') (8) 

R n=l ^ J \ ^ J 

with the interaction tensor 

Wi.iR) = Q{R - a) [fnoiR)5^J + hi2{Rm, - S^i?,)] . (9) 

The functions /no and /112 are certain expansion coefficients of the Mayer function that are 
known analytically. To lowest order in one has /no ~ and fii2{f) ~ {l — e~l^^^')m?/r^ 
so that tfij(R — R') reproduces the dipolar interaction between the points R and R'. This 
expression is minimized with respect to the values of M(R) at all lattice points using the 
simulated annealing algorithm ||69|. Up to = 24 points for each spatial dimension are used 
which amounts to 3 x 24^ = 41472 variables. 

Results for the thermodynamic parameters m* = 1.5, T* = 2.25, and p* = 0.94 are shown 
in Fig. 1^. One of the three perpendicular cubic axes is chosen spontaneously as a vortex axis 
around which the magnetization field circulates. As one can infer from the sections parallel 
to the vortex axis (i.e., the z axis) in Figs. |^(b) and (c) M(r) lies essentially in the vortex 
plane except close to the axis where it "escapes" into the third dimension in order to avoid 
a line singularity. Since the magnetization is approximately parallel at the surface two point 
defects at the centers of the top and bottom face are inevitable. The absolute value M(r), 



which in contrast to the common assumption in micromagnetic calculations |7^ is not fixed 
at a constant value, is strongly reduced within the cores of the point defects whereas it is 
approximately constant throughout the rest of the volume with the exception of thin surface 
layers. At higher temperatures in the vicinity of the phase transition to the isotropic liquid 
a reduction of M(r) occurs also along the whole vortex axis. 

The structure in Fig. ^ can also be considered to consist of four triangular domains with 
four broad 90° domain walls along the diagonals. These walls are clearly induced by the cube 
edges and would probably be absent in a cylindrical sample. We have computed a suitably 
defined wall thickness 6 for different sample sizes L. In Fig. ^ we show the results as 
function of 1/L up to the largest systems with L/a = 12. The various points for the same L 
correspond to different values of A^ and thus give an indication of the error induced by the 
discretization. The data are consistent with the conjecture S/L const for L ^ 00 which 
would mean that the wall thickness is determined by the system size [71]. This contrasts 
with the finite wall thickness in solid ferromagnets that is caused by the existence of easy 
axes induced by the lattice anisotropy. The whole magnetization configuration fulfills, to a 
good approximation, the scaling relation 

M{r/a,L/a) ^ M(°)(r/L) (10) 

when different sizes within the accessible regime are compared. We conjecture that this holds 
also in the thermodynamic limit. 

The average magnetization of the sample (M) = J (i^rM(r) is approximately zero 
although |M(r)| 7^ 0. The positive contribution to (M^) from the "escape" near the vortex 
axis is compensated by negative contributions from the regions close to the edges. However, 
a finite value of (M^) is expected for cylindrical or spherical samples. 
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IV. EXTERNAL FIELD 



An external magnetic field tends to rotate all dipoles into the field direction and thus 
to destroy any inhomogeneity of M(r). Let us first consider sample shapes which induce a 
planar zero-field structure. As we have seen, this is approximately true for the cube, but most 
probably also for a torus [|73[ . We consider the case that the field is applied perpendicularly to 
the magnetization plane. If one assumes that the angle by which the magnetization is rotated 
out of plane is the same throughout the sample one can still treat the inhomogeneous problem 
quasi-analytically |2l|]. Upon increasing the external field the parallel component of M first 
increases linearly and finally saturates while the perpendicular components decrease and 
vanish with a square root singularity at a critical field strength He above which the sample 
is homogeneously magnetized. One can also determine phase diagrams at fixed external 



fields which now do depend on the sample shape via the demagnetization factor D [21 



The results for a cubic or spherical sample [D = 1/3) with m* = 1.5 are shown in Fig. ^. 
The critical line between homogeneous isotropic and inhomogeneous ferromagnetic fluids in 
zero field turns into a critical line between homogeneously and inhomogeneously magnetized 
fluids. With increasing field strength this line together with the tricritical point shifts to 
higher densities and lower temperatures. (Only for a needle-shaped sample {D = 0) this line 
disappears and the tricritical point turns into a second critical point because all phases are 
homogeneous.) The small shift of the liquid-gas critical temperature is proportional to 
for small fields in accordance with the results of perturbation theory |]75[. The influence of 
an external field on the liquid-gas coexistence has also been studied in Refs. [^,0,^]. 

If the external field is not perpendicular to the magnetization plane, for small fields 
the structure is still similar to the zero field result. On the other hand, for example in a 
cube there are six degenerate zero-field solutions. By using the free minimization procedure 
described in the previous section one can address the question which of them is selected by 
a small external field parallel to one of the cube axes. It turns out that a slightly perturbed 
vortex structure with its axis parallel to the field is favored [Q. The further evolution of the 
structure with increasing field strength takes place as described above. Again the azimuthal 
magnetization component around the field direction vanishes at a critical field strength He, 
leaving an almost homogeneously magnetized sample for H > H^. A comparison between 
the structures in zero field and in a finite field below He is given in Figs. |^ and ^. It is 
interesting to consider the corresponding behavior of M(r) for a spherical sample which 
possesses infinitely many equivalent directions. Fredkin and Koehler [74| have performed 
micromagnetic calculations for solid spherical particles. For low lattice anisotropy they find 
again a vortex state, but the vortex axis rotates away from the field direction as H vanishes. 



V. SOLID PHASES 



It is important to obtain an estimate for the position of the freezing transition in the phase 
diagram in order to know whether the formation of the ferromagnetic liquid is preempted 
by freezing. In order to get a reasonable description of the solid phase we have chosen a 
more sophisticated density functional |^. The hard sphere reference fluid is now treated 
within the modified weighted density approximation (MWDA) introduced by Denton and 
Ashcroft which gives very accurate results for the freezing properties of hard spheres. 
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The perturbative contribution is changed in a way that reproduces the results of Curtin and 
Ashcroft in the hmit of a Lennard- Jones fluid but also includes a contribution from the 
Mayer function of the dipolar potential as in Eq. (Q). The number density is parametrized 
as 

where 1 / y/X is the width of the density peaks at the lattice points R. The shape dependent 
contribution to the density functional has been calculated also for a solid structure. For 
cubic lattices it has been shown to be the same as that of a liquid. 

For the considered parameter values [m* < 2) the fee lattice was found to be more 
stable than the bet (body centered tetragonal) structure observed in simulations The 
phase diagrams obtained by this approach (Figs. |^ and H) comprise five different phases, 
among them a ferromagnetic solid (fs) magnetized in the (111) direction, and a plastic solid 
(s) without orientational order. The transition between the liquid phases with and without 
orientational order extends into the solid region. For dipole moments larger than m* ^ 1 
there is indeed a region with a stable ferromagnetic liquid (fl) phase. However, based on 
a different density-functional theory Klapp and Forstmann [|ll]] have claimed recently that 



the latter phase is metastable with respect to the solid for all dipole moments, at variance 
with the simulation results. Thus even the qualitative features of the Stockmayer fluid phase 
diagram are not yet finally settled. 



VI. NONSPHERICAL PARTICLES 



In fluids consisting of nonspherical particles an orient at ionally ordered liquid phase may 
already occur in the absence of dipolar interactions. This nematic phase can be stabilized 
by short-ranged steric interactions between sufficiently elongated or disklike particles. It 
has a higher symmetry than a ferromagnetic liquid as particle orientations parallel and an- 
tiparallel to the preferential direction have identical probabilities. It is interesting to study 
the interplay of these two ordering mechanisms in a system of nonspherical dipolar parti- 
cles. Specifically we have considered hard ellipsoids and hard spherocylinders with a central 
longitudinal dipole moment [^]. In this case, the reference part of the density functional 
coincides with the hard core part which is well approximated by the so-called decoupling 
approximation that amounts to replacing the direct correlation function by that of a hard- 
sphere fiuid scaled by the distance of closest approach for given orientations of the particles 
and their joining vector. The dipolar interaction is again treated perturbatively within the 
low-density approximation that yields an expression bilinear in the densities with the Mayer 
function kernel as in Eq. (|^). 

The properties of this model are determined by the packing fraction rj, the dimension- 
less ratio kBTa^/rn? , where and a\\ are the particle diameters measured perpendicular 
and parallel to the dipole moment, respectively, and the aspect ratio k = (T\\/a_[_ of the 
particles. We find that for elongated particles the isotropic-nematic transition is shifted to 
lower densities by the dipole moment. This trend is agreement with earlier results from the 
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Onsager virial theory [^], integral equation theory ||82|, and a different version of density- 



functional theory |^3[. The two-phase region widens into gas-nematic coexistence at high 
dipole moments or low temperatures (Fig. ||). For ellipsoids of aspect ratio k ^ 2 (Fig. ^) a 
ferromagnetic fluid (F) appears at high densities, as well as a reentrant nematic phase at high 
dipole moments that most probably is an artefact of the theoretical approximations. The 
ferromagnetic phase is reached from the nematic or isotropic states by continuous (dotted 
lines) or weakly first-order (solid lines) transitions, depending on the temperature and par- 
ticle aspect ratio. Similar phase diagrams are found for oblate particles {k < 1, see Fig. ^Ij). 
If one compares fluids with aspect ratios k and 1/k the isotropic-nematic coexistence pack- 
ing fractions are the same for m = (or T —>■ oo) but the ferromagnetic phase becomes 
stable at much higher values of /csTVo/m^ (Vq is the particle volume) for a given packing 
fraction if k < 1, because in this case the dipole moments can approach each other closer in 
the most favorable head-to-tail configuration. For the same reason the ferromagnetic liquid 
should be easier to attain in fluids of oblate than of spherical particles. For dipolar hard 
sphero cylinders one finds an analogous series of phase diagrams as for elongated ellipsoids 



VII. GAS-LIQUID INTERFACE 

Up to this point we were only concerned with bulk properties of dipolar fluids and the 
influence of sample surfaces. If one allows for the presence of interfaces between coexisting 
bulk phases various new phenomena arise. By imposing appropriate boundary condition in 
principle one can examine interfaces between any two phases that can coexist according to 
the phase diagrams shown in the previous sections. These interfaces can be studied with 
the same density-functional theories, although the calculations will be much more involved 
due to the lower symmetry of p{r,u)). Here we restrain ourselves to the simplest case, i.e., 
the interface between the isotropic gas and the isotropic liquid For this system 



the number density depends on the coordinate z normal to the interface and the local 
orientational order near the interface is described by the first nontrivial coefficient a2{z) in 
an expansion in terms of Legendre polynomials: 

p(r, u) = p{z, cos 9) = ^p{z) Q + a2{z)P2{cos 9) + . . }j . (12) 

(The higher coefficients aiy2{z) are very small.) The functions p{z) and a2{z) have been 
calculated numerically by solving the integral equation that results from the stationary 
condition for the density functional. The density profiles are monotonic and for small reduced 
dipole moment m* they depend only weekly on m* for fixed distance r = 1 — T/T^ from 
the critical temperature Tc (see Fig. |1^). Near the interface the particles are preferentially 
aligned parallel to the interface (02 (-z) < 0) on the liquid side and perpendicular to the 
interface (02(2;) > 0) on the vapor side (Fig. p!3| ). The amplitude of 02 (^) increases strongly 
(~ m*^) with the dipole moment. Near the critical point both the density profile 

p*^,^ r^O) = pl + A^im*yF^iz/0 (13) 

and the orientational profile 
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«2(^,r ^ 0) = A^,{7rfy+^^F^,{z/0 (14) 

exhibit a scaling behavior. Here ^ is the bulk correlation length, Fp and F^^ are scaling func- 
tions, and Ap and are nonuniversal amplitudes. Within the density functional approach 
the critical exponents (3 and v of the order parameter and of the correlation length take on 
their mean- field value |. Equation (|1^ has been used to interpret ellipsometric measure- 
ments of the liquid-liquid interface of the critical ionic mixture triethyl-n-hexylammonium 



triethyl-n-hexylboride dissolved in diphenyl ether ||85[| . 

The temperature dependence of the surface tension 7, which can be calculated from these 
profiles, is shown in Fig. lA. The limiting power law 7(r — > 0) = 7oT'^''^ is valid even quite 



far away from T^. The amplitude 70 increases with increasing dipole strength. The values 
of 7 are (for m* < 2) practically unaltered if the orientational profile 02 (-2) is replaced by 
a2{z) = 0. This means that the main influence of the dipolar interaction on the surface 
tension is due to the change of the bulk phase diagram and of the number density profile 
but not due to preferential orientations at the interface. This result has been exploited for 
obtaining estimates of the surface tension of polar fluids from other, more simple analytical 



approaches and simulations 



11 



REFERENCES 



[1] M. Born, Sitz. Phys. Math. 25, 614 (1916); Ann. Phys. 55, 221 (1918). 
[2] D. Wei and G.N. Patey, Phys. Rev. Lett. 68, 2043 (1992). 
[3] D. Wei and G.N. Patey, Phys. Rev. A 46, 7783 (1992). 

[4] J.J. Weis, D. Levesque, and G.J. Zarragoicoechea, Phys. Rev. Lett. 69, 913 (1992). 

[5] J.J. Weis and D. Levesque, Phys. Rev. E 48, 3728 (1993). 

[6] M.J. Stevens and G.S. Grest, Phys. Rev. E 51, 5962 (1995). 

[7] M.J. Stevens and G.S. Grest, Phys. Rev. E 51, 5976 (1995). 

[8] K. Sano and M. Doi, J. Phys. Soc. Jpn. 52, 2810 (1983). 

[9] H. Zhang and M. Widom, Phys. Rev. E 49, R3591 (1994). 
[10] D. Wei, G.N. Patey, and A. Percra, Phys. Rev. E 47, 506 (1993). 
[11] S. Klapp and F. Forstmann, Europhys. Lett. 38, 663 (1997). 
[12] H. Zhong and R.G. Petschek, Phys. Rev. E 51, 2263 (1995). 
[13] J.S. H0ye and G. Stell, MoL Phys. 86, 707 (1995). 
[14] S. Klapp and F. Forstmann, J. Chem. Phys. 106, 9742 (1997). 
[15] M.E. van Leeuwen, B. Smit, and E.M. Hendriks, MoL Phys. 78, 271 (1993). 
[16] G. Kronome, J. Liszi, and L Szalai, J. Chem. Soc. Faraday Trans. 93, 3053 (1997); and 
references therein. 

[17] L Szalai, G. Kronome, and T. Lukacs, J. Chem. Soc. Faraday Trans. 93, 3737 (1997). 
[18] C. Kricbel and J. Winkelmann, Mol. Phys. 90, 297 (1997); and references therein. 
[19] B. Groh and S. Dietrich, Phys. Rev. Lett. 72, 2422 (1994); ibid 74, 2617 (1995); Phys. 

Rev. E 50, 3814 (1994). 
[20] M.A. Osipov, P.LC. Teixeira, and M.M. Telo da Gama, J. Phys. A. 30, 1953 (1997). 
[21] B. Groh and S. Dietrich, Phys. Rev. E 53, 2509 (1996). 
[22] R.B. Griffiths, Phys. Rev. 176, 655 (1968). 

[23] J.M. Tavares, M.M. Telo da Gama, P.I.C. Teixeira, J.J. Weis, and M.J.P. Nijmeijer, 

Phys. Rev. E 52, 1915 (1995). 
[24] M.J.P. Nijmeijer and J.J. Weis, Ann. Rev. Comp. Phys. IV, 1 (1996). 
[25] A. Oukouiss and M. Bans, Phys. Rev. E 55, 7242 (1997). 
[26] M.J.P. Nijmeijer and J.J. Weis, Phys. Rev. Lett. 75, 2887 (1995). 
[27] M.J.P. Nijmeijer and J.J. Weis, Phys. Rev. E 53, 591 (1996). 
[28] M.J.P. Nijmeijer, A. Parola, and L. Reatto, Phys. Rev. E 57, 465 (1998). 
[29] M.A. Zaluska-Kotur and L.A. Turski, Phys. Rev. A 41, 3066 (1990). 
[30] R.E. Roscnswcig, Ferrohydrodynamics (Cambridge University Press. Cambridge, 1985); 

E. Blums, A. Ccbers, and M.M. Maiorov, Magnetic Fluids (de Gruyter, Berlin, 1997). 
[31] Proceedings of the 5th International Conference on Electro-Rheological Fluids, 

Magneto-Rheological Suspensions and Associated Technology, Int. J. Mod. Phys. B 

10, Numbers 23 & 24 (1996). 
[32] D. Platzck, C. NotthofT, D.M. Herlach, G. Jacobs, D. Herlach, and K. Maier, Appl. 

Phys. Lett. 65, 1723 (1994). 
[33] J. Reske, D.M. Herlach, F. Keuser, K. Maier, and D. Platzek, Phys. Rev. Lett. 75, 737 

(1995). 

[34] T. Albrecht, C. Biihrer, M. Fahnle, K. Maier, D. Platzek, and J. Reske, Appl. Phys. A 

65, 215 (1997). 

[35] J.J. Weis and D. Levesque, Phys. Rev. Lett. 71, 2729 (1993). 



12 



D. Levesque and J.J. Weis, Phys. Rev. E 49, 5131 (1994). 

M.E. van Leeuwen and B. Smit, Phys. Rev. Lett. 71, 3991 (1993). 

J.M. Tavares, M.M. Telo da Gama, and M.A. Osipov, Phys. Rev. E 56, R6252 (1997). 

S.C. McGrother and G. Jackson, Phys. Rev. Lett. 76, 4183 (1996). 

R.P. Sear, Phys. Rev. Lett. 76, 2310 (1996). 

R. van Roij, Phys. Rev. Lett. 76, 3348 (1996). 

M.A. Osipov, P.I.C. Teixeira, and M.M. Tclo da Gama, Phys. Rev. E 54, 2597 (1996). 

T. Chen, R.N. Zitter, and R. Tao, Phys. Rev. Lett. 68, 2555 (1992). 

W.E. McMuUen and D.W. Oxtoby, J. Ghem. Phys. 86, 4146 (1987). 

S.J. Smithhne, S.W. Rick, and A.D.J. Haymet, J. Ghem. Phys. 88, 2004 (1988). 

W.E. McMullen and D.W. Oxtoby, J. Ghem. Phys. 88, 1476 (1988). 

B. Groh and S. Dietrich, Phys. Rev. E 54, 1687 (1996). 

B.E. Vugmeister and M.D. Ghnchuk, Rev. Mod. Phys. 62, 993 (1990). 

H. Zhang and M. Widom, J. Magn. Magn. Mat. 122, 119 (1993); Phys. Rev. B 51, 8951 

(1995). 

G. Ayton, M.J.P. Gingras, and G.N. Patcy, Phys. Rev. Lett. 75, 2360 (1995). 

G. Ayton and G.N. Patey, Phys. Rev. Lett. 76, 239 (1996). 

J. Eggebrecht, K.E. Gubbins, and S.M. Thompson, J. Ghem. Phys. 86, 2286 (1987). 

M. Kasch and F. Forstmann, J. Ghem. Phys. 99, 3037 (1993). 

P.L Teixeira and M.M. da Gama, J. Phys. Gondens. Matter 3, 111 (1991). 

M. Matsumoto, H. Mizukuchi, and Y. Kataoka, J. Ghem. Phys. 98, 1473 (1993). 

N.A.M. Bessehng and J. Lyklcma, J. Phys. Ghem. 98, 11610 (1994). 

M.J. Booth, D.-M. Duh, and A.D.J. Haymet, J. Ghem. Phys. 101, 7925 (1994). 

B. Yang, D.E. SuUivan, and G.G. Gray, J. Phys. Gondens. Matter 6, 4823 (1994). 

V.P. Sokhan and D.J. Tildesley, MoL Phys. 92, 625 (1997); and references therein. 

V. Talanquer and D.W. Oxtoby, J. Ghem. Phys. 99, 4670 (1993). 

V. Talanquer and D.W. Oxtoby, J. Ghem. Phys. 103, 3686 (1995). 

H. B. Lavender, K.A. Iyer, and S.J. Singer, J. Ghem. Phys. 101, 7856 (1994). 
D. Lu and S.J. Singer, J. Ghem. Phys. 103, 1913 (1995). 

D.J. Henderson and W. Schmickler, J. Ghem. Soc. Faraday Trans. 92, 3839 (1992). 
G. Kricbel and J. Winkelmann, MoL Phys. 90, 297 (1997); and references therein. 
P. Frodl and S. Dietrich, Phys. Rev. A 45, 7330 (1992); Phys. Rev. E 48, 3203 (1993). 
J. A. Barker and D. Henderson, J. Ghem. Phys. 47, 4714 (1967). 

B. Groh and S. Dietrich, Phys. Rev. Lett. 79, 749 (1997); Phys. Rev. E, 57, 4535 (1998). 
W.H. Press, B.P. Flannery, S.A. Teukolsky, and W.T. Vetterhng, Numerical Recipes 
(Cambridge University Press, Cambridge, 1989). 

W.F. Brown jr., Micromagnetics (Krieger, Huntington, 1978). 
P.G. dc Gcnnes and P.A. Pincus, Sohd State Communications 7, 339 (1969). 
A. Aharoni and J.P. Jakubovics, J. Magn. Magn. Mat. 83, 451 (1990). 
P.J. Wojtowicz and M. Rayl, Phys. Rev. Lett. 20, 1489 (1968). 

D.R. Fredkin and T.R. Koehler, IEEE Trans. Magn. MAG-24, 2362 (1988); J. AppL 
Phys. 67, 5544 (1990). 

D. Boda, I. Szalai, and J. Liszi, J. Ghem. Soc. Faraday Trans. 91, 889 (1995). 

C. E. Woodward and S. Nordholm, J. Phys. Ghem. 92, 501 (1988). 

D. Boda, J. Winkelmann, J. Liszi, and I. Szalai, Mol. Phys. 87, 601 (1996). 



13 



[78] A.R. Denton and N.W. Ashcroft, Phys. Rev. A 39, 4701 (1989). 

[79] W.A. Curtin and N.W. Ashcroft, Phys. Rev. Lett. 56, 2775 (1986). 

[80] B. Groh and S. Dietrich, Phys. Rev. E 55, 2892 (1997). 

[81] A.G. Vanakaras and D.J. Photinos, MoL Phys. 85, 1089 (1995). 

[82] A. Perera and G.N. Patey, J. Chem. Phys. 91, 3045 (1989). 

[83] C. Vega and S. Lago, J. Chem. Phys. 100, 6727 (1994). 

[84] P. Frodl and S. Dietrich, Phys. Rev. E 48, 3741 (1993). 

[85] C.L. Caylor, B.M. Law, P. Senanayake, V.L. Kuzmin, V.P. Romanov, and S. Wiegand, 

Phys. Rev. E 56, 4441 (1997). 
[86] S. Abbas, P. Ahlstromega, and S. Nordholm, Langmuir 14, 396 (1998). 



14 



FIGURES 



FIG. 1. Phase diagrams of the Stockmayer fluid for three different values of the reduced dipole 
moment m*, including the ferromagnetic liquid phase. The dashed line pfdT) denotes a line of 
critical points, Tc is the gas-liquid critical point, and the shaded areas represent two-phase regions. 
In part (a) the critical line ends at the coexistence line in a critical end point Tcep, while in (b) 
and (c) a tricritical point Tt arises, below which the isotropic-ferromagnetic transition becomes 
first order. The dotted lines represent the metastable parts of the coexistence curve of the isotropic 
phases. 

FIG. 2. Magnetization structure of a ferromagnetic liquid in a cubic volume with edge length 
L/a = 7.2, 20^ mesh points, and T* = 2.25. The three parts of the figure represent sections (a) for 
fixed z = 0.025L, (b) fixed y = — 0.075L, and (c) fixed y = — 0.225L; the origin is at the center of 
the cube. The arrows represent projections of the local magnetization at their midpoint onto the 
section plane. The z axis is the vortex axis. Part (c) demonstrates that in the vicinity of this axis 
the magnetization "escapes" from the vortex plane. 

FIG. 3. The domain wall thickness 5 in a cubic sample measured within the plane z = at 
distance r = L/2^/^ from the center, i.e., at half distance between the center and the edges. The 
different results for the same value L* = L/a correspond to different numbers of mesh points 
used for the minimization. The extrapolation to 1/L = suggests a finite value of 6/L in the 
thermodynamic limit. 

FIG. 4. Phase diagrams for m* = 1.5 and various fixed values of the reduced external field 
strength H* = H^Ja^Je for a sample with aspect ratio k = \ (i.e., demagnetization factor D = 1/3). 
For nonzero fields the inhomogeneous domain structure is partly retained at high densities and low 
temperatures. The magnetic field shifts the transition line of critical points (dotted curve) and 
the tricritical point (open circle) towards higher densities and lower temperatures, whereas the 
gas-liquid critical temperature (full circle) is slightly increased. 

FIG. 5. Sections of the magnetization structure orthogonal to the field direction for L/a = 9.6 
and N = 16 at z/a = —0.9 (a) without and (b) with an external field H* = H-\fa^e = 1 applied 
in the z direction. The scale factor which determines the lengths of the arrows is the same in both 
parts of this figure and also in Fig. |^ below. The transversal components of the magnetization are 
reduced by the field but the overall feature of the inhomogeneous structure is preserved. 

FIG. 6. Vertical sections of the same magnetization structures as in Fig. |5| perpendicular to the 
y axis and thus parallel to the field direction for y/a = —2.1. At this distance of the plane from 
the center the z component of the magnetization in zero field is only small, while the external field 
induces a large z component everywhere. 
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FIG. 7. Phase diagram for m* = 1 including the sohd phase as obtained from the generahzed 
density- functional theory. It encompasses five phases: the isotropic gas (g) and liquid (1), the 
ferromagnetic liquid (fl), the ferromagnetic solid (fs), and the solid without orientational order 
(s). Two-phase regions are shaded, the dashed lines indicate continuous phase transitions, and the 
dotted line denotes the triple point. Critical end points are denoted by squares. 

FIG. 8. The same as Fig. |^ but for m* = 1.35. In addition to the features of Fig. ^ a tricritical 
point (open circle) and two triple points (dotted lines) occur. 

FIG. 9. Phase diagram of hard ellipsoids of revolution with a longitudinal dipole moment and 
aspect ratio k = 3. rj is the packing fraction. At high temperatures or small dipole moments the 
transition between the isotropic (I) and the nematic (N) phases is weakly first order, but it widens 
into a gas-nematic coexistence at low temperatures. Possible smectic and solid phases cannot be 
described within the present theory. 

FIG. 10. The same as Fig. ^ but for aspect ratio k = 2. Besides the isotropic (I) and nematic (N) 
phases a ferromagnetically ordered liquid (F) occurs in the intermediate temperature range. The 
dotted and solid lines denote second- and first-order phase transitions, respectively. The density 
gap of the I-N transition at 77 ~ 0.68 cannot be resolved on the present scale. The inset shows that 
in a narrow temperature range the high temperature continuous N-F transition is turned into a 
weakly first-order transition generating a tricritical point and an I-N-F triple point. The reentrant 
nematic phase at low temperatures is probably an artefact of the present theory. 

FIG. 11. Phase diagram for oblate ellipsoids with aspect ratio k = 1/3. The I-N coexistence 
packing fractions for T ^ 00 are the same as those for k = 3, but in contrast to Fig. ^ a ferromag- 
netic phase is stable at high densities. 

FIG. 12. Density profiles for the gas-liquid interface of a Stockmayer fluid for various dipole 
moments at a fixed reduced temperature r = 1 — T/Tc = 0.2. 

FIG. 13. Orientational order profiles 02(2;) at the gas-liquid interface for the same parameters 
as in Fig. |l2l Note that 0:2 is scaled by m*'^. The dipoles are aligned preferentially parallel to the 



interface on the liquid side and perpendicular to the interface on the gas side. For m* — > the 
profiles approach a limiting form a2{z/a,m*) = m*^a^2\^/'^)- 

FIG. 14. Temperature dependence of the liquid-gas surface tension 7^*^ = ji^gcr'^/e for various 
dipole moments. The double logarithmic plot in the inset demonstrates that ji^g ~ r^/^ for r — > 0. 
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